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Q^ \ Abstract 

, The solution of QCD equations for generating functions of parton 

Q I muhiphcity distributions reveals new peculiar features of cumulant 

' moments oscillating as functions of their rank. It happens that ex- 

perimental data on hadron multiplicity distributions in e~^e~ ,hh, AA 
collisions possess the similar features. However, the "more regular" 
models like Xcp^ behave in a different way. Evolution of the moments 
^ . at smaller phase space bins and zeros of the truncated generating 

^T) ' functions are briefly discussed. 
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1 Introduction and main results 

Since the space-time (space in proceedings and time at oral presentation) is 
' very limited I present from the beginning both short review of the problem 

and the main results obtained, leaving their justification and brief discussion 
for the next section. Those interested in more detailed description should 
^ ■ use the list of references for further reading (in particular, the review paper 

^ i I would like to stress that QCD predicts the distributions of partons 

^ I (quarks and gluons) while in experiment one gets the distributions of final 

hadrons. Therefore no quantitative comparison has been attempted. How- 
ever, the qualitative features of both distributions are so spectacular and 
remind each other that one is tempted to confirm once again that QCD is a 
powerful tool for predicting new features of hadron distributions as well. 

For a long time, the phenomenological approach dominated in descrip- 
tion of multiplicity distributions in multiparticle production The very 
first attempts to apply QCD formalism to the problem failed because in 
the simplest double-logarithmic approximation it predicts an extremely wide 
shape of the distribution that contradicts to experimental data. Only 
recently it became possible to get exact solutions of QCD equations for 
the generating functions of multiplicity distributions which revealed much 



narrower shapes and such a novel feature of cumulant moments as their os- 
cillations at higher ranks 0,0. The similar oscillations have been found in 
experiment for the moments of hadron distributions [0]. Their pattern differs 
drastically from those of the popular phenomenological distributions |l| and 
of the "non-singular" A0^-model 

These findings have several important implications They show that: 
l.the QCD distribution belongs to the class of non-infinitely- divisible ones. 
Two corollaries of this statement follow immediately: 

a. The Poissonian cluster models (e.g., the multiperipheral cluster model) 
are ruled out by QCD. 

b. The negative binomial distribution (so popular nowadays in phenomeno- 
logical fits) can not be valid asymptotically. 



2. the new expansion parameter appears in description of multiparticle pro- 
cesses. 

Since this parameter becomes large when large number of particles are in- 
volved, it asks for the search of some collective effects and of more convenient 
basis than the common particle number representation. 

QCD is also successful in qualitative description of evolution of multiplic- 
ity distributions with decreasing phase space bins which gives rise to notions 
of intermittency and fractality [p|,|IO[ 



nfl. However, there are some new 



problems with locations of the minimum of cumulants at small bins [12|,[13 



The experimentally defined truncated generating functions possess an 
intriguing pattern of zeros in the complex plane of an auxiliary variable [13], 
|l^,[jl5l. It recalls the pattern of Lee- Yang zeros of the grand canonical 



partition function in the complex fugacity plane related to phase transition 
and asks for some collective effects to be searched for UlSj , ||T9[| . 



2 Technicalities 

Let us define the multiplicity distribution 

oo 

n=0 



where (T„ is the cross section of n-particle production processes, and the 
generating function 

oo 

G{z) = Y.Pn{l + zr. (2) 

n=0 

The (normahzed) factorial and cumulant moments of the P„ distribution are 

J, _ T.nPnn{n-l)...{n-q + l) 1 d'iG{z) 

_ 1 dnnGjz) 
^''W dz^ '^=°' 
where (n) = J2n is the average multiphcity. They describe full and 
genuine g-particle correlations, correspondingly. Let us point out here that 
the moments are defined by the derivatives at the origin and are very sensitive 
to any nearby singularity of the generating function. 

First, let us consider QCD without quarks, i.e. gluodynamics. The gen- 
erating function of the gluon multiplicity distribution in the full phase-space 
volume satisfies the equation 

= dxK{x)^l[G{z, y + In y + ln(l - x)) - G{z, Y)] . (5) 

Here Y = \n{p6/Qo),p is the initial momentum, 6 is the angular width of the 
gluon jet considered, p9 = Q where Q is the jet virtuality, Qo =const, 

2 6as{Q) , . 

To = — - — , (6) 

TT 

as is the running coupling constant, and the kernel of the equation is 

K{x) = - -{l-x)[2-x{l-x)]. (7) 

X 

It is the non-linear integro-differential equation with shifted arguments in the 
non-linear part which take into account the conservation laws. , and with 
the initial condition 

G{z,Y = Q) = l + z, (8) 

and the normalization 

G{z = 0,Y) = l. (9) 



The condition (|^) normalizes the total probability to 1, and the condition 
(H) declares that there is a single particle at the very initial stage. 

After Taylor series expansion at large enough Y and differentiation in eq. 
(H), one gets the differential equation 

(InG(F))" = 7o'[G(r) - 1 - 2hG'{Y) + hG"iY)], (10) 

where hi = 11/24; /i2 = (67 — 67r^)/36 ~ 0.216, and higher order terms have 
been omitted. 

Leaving two terms on the right-hand side, one gets the well-known equa- 
tion of the double-logarithmic approximation which takes into account the 
most singular components. The next term, with hi, corresponds to the mod- 
ified leading-logarithm approximation, and the term with /i2 deals with next- 
to-leading corrections. 

The straightforward solution of this equation looks very problematic. 
However, it is very simple for the moments of the distribution because G{z) 
and lnG{z) are the generating functions of Fq and Kg, correspondingly, ac- 
cording to (^, (^. Using this fact, one gets the solution which looks like 

rr 7g[l-2/iig7 + /^2(gV + g70] 

where the anomalous dimension 7 is related to 70 by 

7 = 70 - lha', + 1(4/^2 - hl)^, + 0(7o^). (12) 

The formula (pH]) shows how the ratio Hg behaves in different approxi- 
mations. In double-log approximation, where hi = h2 = 0, it monotonously 
decreases as that corresponds to the negative binomial law with its pa- 
rameter k = 2 i.e. to very wide distribution. Let us note that it owes to the 
singular part of the kernel and is absent on more regular theories like X(j)^. 
In modified-log approximation {h2 = 0) it acquires a negative minimum at 

Qmin = r— + l+ 0(70) ~ 5 (13) 

/ii7o 2 

and approaches the abscissa axis from below asymptotically at large ranks 
q. In the next approximation given by ([TTl) it preserves the minimum loca- 
tion but approaches a positive constant crossing the abscissa axis. In ever 



higher orders it reveals the quasi-oscillatory behavior about this axis. This 
prediction of the minimum at g ~ 5 and subsequent specific oscillations is 
the main theoretical outcome. 

It is interesting to note that the equation (|) can be solved exactly in 
the case of fixed coupling constant . All the above qualitative features are 
noticeable here as well. 

While the above results are valid for gluon distributions in gluon jets (and 
pertain to QCD with quarks taken into account the similar qualitative 
features characterize the multiplicity distributions of hadrons in high energy 
reactions initiated by various particles. The numerous demonstration of it 
can be found in the review paper 

The multiplicity distributions can be measured not only in the total phase 
space (as has been discussed above for very large phase-space volumes) but in 
any part of it. For the homogeneous distribution of particles within the vol- 
ume, the average multiplicity is proportional to the volume and decreases for 
small volumes but the fiuctuations increase. The most interesting problem 
here is the law governing the growth of fiuctuations and its possible depar- 
ture from a purely statistical behavior related to the decrease of the average 
multiplicity. Such a variation has to be connected with the dynamics of the 
interactions. In particular, it has been proposed to look for the power-law 
behavior of the factorial moments for small rapidity intervals 5y 



inspired by the idea of intermittency in turbulence. In the case of statistical 
fiuctuations with purely Poissonian behavior, the intermittency indices 0(g) 
are identically equal to zero. 

Experimental data on various processes in a wide energy range support 
this idea , and QCD provides a good basis for its explanation as a result of 
parton showers. 

At moderately small rapidity windows, one can get in the double-log 
approximation the power-law behavior with 



F.^idyr^^"^ (0(g) >0) (51/ -0), 



(14) 




(15) 



The running property of QCD coupling constant is not important in that 
region. This property becomes noticeable at ever smaller windows, where 



(e.g., at q=2) \n6yo/Sy > a^^, and leads to smaller numerical values of 0(g) 
compared to (pIS]). The general trends in this region decline somewhat from 
the simple power law (|T3|) due to logarithmic corrections. Qualitatively, these 
predictions correspond to experimental findings at relatively small ranks q 
where the steep increase in the region of Sy > 1 on the log-log plot of the de- 
pendence (0) is replaced by slower one at smaller intervals 6y. The transition 
point between the two regimes depends on the rank in qualitative agreement 
with QCD predictions also. Namely, the transition happens at smaller bins 
for higher ranks. These findings can be interpreted as an indication on frac- 
tal structure of particle distributions within the available phase space. When 
interpreted in terms of fiuctuations, they show that the fiuctuations become 
stronger in small phase-space regions in a definite power-like manner and, 
surely, exceed trivial statistical fiuctuations. 

Let us turn now to the g-behavior of moments at small bins. The phe- 
nomenon of the oscillations of cumulants discussed above reveals itself here 
as well if one goes beyond the double-log approximation of (|15|). In terms of 
factorial moments, it means the non-monotonous behavior of the intermit- 
tency indices as functions of q. (Compare it to the steady increase with q at 



g > 1 given by (|T5|).) It gives rise to the negative values of Kg and Hq. 

The fate of the first minimum can be easily guessed from the formula 
(|T3|). For large enough virtualities (i.e. small 70), the minimum location 
moves to higher values of rank g for jets with larger virtuality Q since the 
QCD coupling constant is running as ln~^ Q. Therefore, the predicted shift 
of the minimum is 

qmin oc In^/^ Q. (16) 

It follows that qmin moves to higher ranks at higher energies because more 
massive jets become available. Another corollary is that it should shift to 
smaller values of g for smaller bins at fixed energy. 

While former statement finds some support in experiment, the second 
one does not look to be true. On the contrary, the minimum appears at 
higher ranks for smaller bins. There is no solution of this problem yet but it 
should be ascribed to the higher-order effects. Actually, one can guess that 
the higher order terms shown as 0(70) in ( |l3l) become so important at small 
bins that they overpower the weak Q-dependence of 'Jq^ in the first term 
of ([T3|) . It is important to stress here that at large rapidity intervals the 
modified leading-log term with h2 does not infiuence the value of qmin, and 



only increases the value of Hq by /i27o- Thus, the next- to-leading corrections 
should be in charge of the additional shift of gmm, and, therefore, small bins 
help us look into higher orders of QCD. 

There is another fascinating feature of multiplicity distributions - it hap- 
pens that zeros of the truncated (sum in || runs up to n = N^ax) generating 
function form a spectacular pattern in the complex plane of the variable z. 
Namely, they seem to lie close to a single circle. At enlarged values of N^ax 
they move closer to the real axis pinching it at some positive value of z. 

No QCD interpretation of the fact exists because it is hard to exploit 
the finite cut-off in analytic calculations. The interest to it stems from the 
analogy to the locations of zeros of the grand canonical partition function as 
described by Lee and Yang who related them to possible phase transitions 
in statistical mechanics. In that variable plays the role of fugacity, 

and pinching of the real axis implies existence of two phases in the system 
considered. 

In particle physics, it shows up the location of the singularity of the 
generating function i.e. the number of zeros of truncated generating functions 
increases and they tend to move to the singularity point when N^ax — *■ 
oo. Since it happens to lie close to the origin, it drastically influences the 
behavior of moments (see (^), (P), and, therefore, determines the shape of 
the distribution. The study of the singularities is at the very early stage now, 
and one can only say that the singularity is positioned closer to the origin 
in nucleus-nucleus collisions and it is farthest in e~^e~ that appeals to our 
intuitive guess. 

To conclude, I would like to stress that, once again, QCD demonstrates 
its power in predicting new features of particle distributions when dealing 
with parton distributions. 
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